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This paper discusses the verification of the Constraint Force Equation (CFE) methodol- 
ogy and its implementation in the Program to Optimize Simulated Trajectories II (POST2) 
for multibody separation problems using three specially designed test cases. The first test 
case involves two rigid bodies connected by a fixed joint; the second case involves two rigid 
bodies connected with a universal joint; and the third test case is that of Mach 7 separation 
of the Hyper-X vehicle. For the first two cases, the POST2/CFE solutions compared well 
with those obtained using industry standard benchmark codes, namely AUTOLEV and 
ADAMS. For the Hyper-X case, the POST2/CFE solutions were in reasonable agreement 
with the flight test data. The CFE implementation in POST2 facilitates the analysis and 
simulation of stage separation as an integral part of POST2 for seamless end-to-end simula- 
tions of launch vehicle trajectories. 
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Nomenclature 

Rigid Body A, and B 

Joint location in Body A, and B 

Components of sensed mass center acceleration in local body frame 

Angle-of-attack 

Angle-of-sideslip 

x, y and z components of joint translational displacement in local body frame 

Inertia tensor about mass center for Body A, and B 

Joint constraint force vector for Body A, and B 

Constraint force components in local body frame 

External force vector acting on Body A, and B 

External force components in local body frame 

Mass of Body A, and B 

Angular velocity vector of Body A, and B relative to an inertial reference frame 
Inertial position vector of Body A, and B 

Vector from Body A mass center to A , and Body B mass center to B 
Joint constraint torque vector for Body A, and B 
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Constraint torque components in local body frame 
External torque vector acting on Body A, and B 

Inertial components of mass center velocity relative to an inertial reference frame 
Body components of inertial angular velocity. 

Inertial position vector to mass center of Body A, and B 


I. Introduction 

Analyzing the dynamic separation of multiple bodies within the atmosphere is complex and challenging. One 
problem that has received significant attention in the literature is that of store separation from aircraft. 1 A similar 
example is the separation of the X-15 research vehicle from the B-52 carrier aircraft. 2 In both of these cases, the 
store and the X-15 vehicle are much smaller in size than the parent vehicle. The other class of stage separation prob- 
lem involves separation of two vehicles of comparable sizes, as in the case of multi-stage reusable launch vehicles 
where the integrity of each stage is important after separation. 

NASA studies on stage separation of multi-stage reusable launch vehicles date back to the early 1960s. 3 " 7 These 
studies addressed the problem of separation of generic two-stage reusable launch vehicles. More recently, Naftel et 
al. considered staging of dual wing-body vehicles. 8 " 10 NASA’s interest in stage separation research was renewed in 
early 2000 when it was realized that the technologies needed for the development of a next generation, reusable sin- 
gle-stage-to-orbit vehicle were not yet available and the focus shifted to multi-stage launch vehicles. Accordingly, 
NASA initiated a comprehensive stage separation tool development activity that included wind tunnel testing as well 
as development and validation of CFD and engineering level simulation tools. 11 As a part of this activity, a stage 
separation analysis and simulation tool was developed called ConSep (short for Conceptual Separation), which is a 
front end to the commercial multibody dynamic analysis software ADAMS®. 12 References 13-15 discuss ConSep, 
and its application to two-body and three-body separation problems. One disadvantage of ConSep is that it cannot be 
easily integrated into standard trajectory simulation software for performing efficient and seamless end-to-end simu- 
lations of launch vehicle trajectories because it is tied to ADAMS. To address this problem, the authors developed 
the CFE methodology and implemented it into the Program to Simulate Optimized Trajectories II (POST2). 16 

The objective of this paper is to test and verify the CFE methodology, presented by the authors in Ref. 16, as a 
means in which any standard trajectory simulation program capable of simulating multiple unconnected vehicles can 
be modified, with a minimal set of changes, to solve constrained motion problems; that is, problems where vehicles 
are connected by simple joints. Although the CFE methodology is generic in nature and can be implemented in any 
conventional trajectory simulation program, this paper discusses its implementation in the industry standard trajec- 
tory simulation and optimization software package, POST2. 17,18 

By itself, POST2 does not have the capability to model internal joint forces and moments prior to separation 
when the bodies are still connected. The CFE implementation in POST2 provides a framework for computing the 
internal constraint forces and moments acting at joints connecting multiple vehicles and applies them as external 
forces/moments on each body, along with the usual external forces and moments due to gravity, aerodynamics, and 
propulsion that are computed by POST2. Thus, the CFE methodology simply augments the external vehicle external 
loads and does not require modification of the POST2 equations of motion. It provides the missing link to accurately 
model the dynamics of generic stage separation problems. 

In this paper the CFE methodology is briefly reviewed, and some of the steps taken to verify and validate its re- 
sults are highlighted. In particular, its application to three test cases of increasing complexity is described. The first 
test case simulates the motion of two rigid bodies connected by a fixed joint with no external forces. The advantage 
of studying this basic problem is that numerical results can be compared with what is to be expected on the basis of 
first principles. The next test case involves comparisons between results obtained with CFE and ADAMS for two 
rigid bodies connected by a universal joint. Finally, the paper concludes with a comparison of a POST2/CFE simula- 
tion of the Mach 7 separation of the X-43A vehicle to previously published POST2 results (generated using an engi- 
neering model) and the flight data. These three test cases were part of a larger, comprehensive series of checks that 
were performed on the CFE technique that provide confidence in its ability to solve constrained motion problems in 
conventional trajectory programs like POST2. For the first two cases, the CFE simulations are in excellent agree- 
ment with simulations using benchmark, industry standard software. For the third test case, CFE performs as well as 
the engineering model but is easier to implement, and the results compare reasonably well with the flight data. 
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II. Constraint Force Equation Methodology 

To illustrate the basic con- 
cept of the CFE methodology, 
consider the motion of two rigid 
bodies connected by a single 
joint as shown in Fig. 1. The 
external forces and moments 
that act on each body are shown 
in Fig. 1(a). These external 
forces ( F] ext) and F^ EXT) ) and 
external moments or torques 
about mass centers ( T^ ext) and 
T] EXT) ) are the resultants of 
gravity, aerodynamic, and pro- 
pulsive forces on each vehicle 
and are readily computed by 
POST2. Figure 1(b) shows the 
internal constraint forces and 
moments (F (CON) and T (CON) ) 
computed by CFE that act on 
each vehicle at the joint loca- 
tion. These internal forces and 
moments constrain the way in which one vehicle can move relative to the other, and are dependent upon the external 
forces acting on each vehicle as well as the type of joint. For joints that permit only relative rotation between two 
bodies, the forces and moments on one body have magnitudes that are equal and directions that are opposite to those 
acting on the other body, as shown in Fig. 1(b). Finally, Fig. 1(c) illustrates the way in which CFE is implemented in 
POST2. At each integration time step the current version of POST2 computes the typical external forces and mo- 
ments acting on each vehicle. This information, along with specific geometric information about the joint, is pro- 
vided to the CFE routine, which works in parallel to compute the internal constraint forces and moments required to 
satisfy the constraints imposed by the joint. Next, these internal joint loads are applied to each vehicle as additional 
external forces and moments, and the POST2 solution is propagated in the usual manner to the next time step. Thus, 
the net external forces and moments on each vehicle are the sum of the usual external forces and moments and the 
joint loads, applied to each vehicle as additional external forces and moments. Consequently, the CFE joint model 
simply augments the vehicle external loads and does not require modification to the POST2 equations of motion. 

The equations of constrained motion of two rigid bodies A and B connected by a single joint are as follows. For 
body A, 



Figure 1. Decomposition of external forces and internal constraint forces to 
free body diagram solved by CFE method. 
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where p A is the position vector from the mass center of A to the point A of A at which the constraint force is applied. 
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Note that in the case of unconstrained motion, the constraint forces and moments vanish; Eqs. (l)-(4) reduce to 
the standard equations of motion for two independent rigid bodies. Each mass center acceleration, rigid body angu- 
lar acceleration, constraint force, and constraint torque are regarded as unknowns; therefore, Eqs. (l)-(4) as they are 
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written represent 12 scalar equations in 24 unknowns. Hence, twelve additional scalar equations are required. Of 
these, six are obtained from two vector relationships that are consequences of the law of action and reaction: 

F* con) + Ff 0N) = 0 (5) 

T (CON) + t (con) + (rfi _ x f <con) = Q (6) 

where r A is the position vector from a point fixed in an inertial reference frame to A . Because x A is the position vec- 
tor from a point fixed in an inertial reference frame to the mass center of A, we can write r A = x A + p A . It should be 
kept in mind that the vector p A is fixed in A. The position vector r B is defined similarly. Note that the cross product 
term involving (r B - r A ) in Eq. (6) is required for joints that allow translation (see Ref. 16); the equation gives a sim- 
ple, expected result when no translation is permitted (r B = r^). 

Equations (l)-(6) are applicable to each of the three test cases discussed in this paper. The six remaining equa- 
tions fall into one of two categories. First, equations of constraint are formed that describe the ways in which a par- 
ticular joint restricts relative translation or relative rotation. Second, by considering the translation or rotation that is 
permitted by the joint, equations are written to account for the absence of constraint force or torque in a particular 
direction in view of the ideal (perfectly smooth) nature of the joint. 

When relative translation is constrained, the distance between two points must remain fixed in a particular direc- 
tion. That is, 


=0 (7) 

where e A is a unit vector fixed in body A in the direction that the translation motion is not permitted by the joint. One 
equation having the form of Eq. (7) is needed to account for each direction that translation is constrained. For exam- 
ple, a fixed joint constrains translation in three orthogonal directions; therefore, three equations in the form of Eq. 
(7) are required. In each such equation, the role of e A is played by one of the three mutually orthogonal unit vectors 
fixed in A. For a prismatic (sliding) joint that permits translation in only one direction (and restricts translation in 
two perpendicular directions), two equations having the form of Eq. (7) are required. 

A constraint on relative rotation can be viewed as a requirement that two unit vectors must remain perpendicular; 
each unit vector is fixed in one of the bodies and is normal to the axis about which rotation would take place if the 
constraint were not present. The constraint is expressed by setting the scalar product of the two unit vectors equal to 
zero. 


e 5 -e^=0 (8) 

where the two unit vectors e A and e B (one fixed in each body) are chosen to be perpendicular to each other through- 
out the constrained motion. One equation in the form of Eq. (8) is required for each direction about which rotation is 
constrained. 

The constraint equations (7) and (8) are differentiated twice with respect to time so that the resulting equations 
involve the unknown linear and angular accelerations of bodies A and B , 

■«.=»] (9) 

fr[V‘,=o] (10) 

The result of differentiating Eq. (7) once is 

(r s -^)-e J +(r s -r^)-((i0^xe^) = 0 (11) 
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where r A and r B are, respectively, the time derivatives in an inertial reference frame of r A and r B . Differentiating 
once more, one obtains 

(x s + (i) B x - x^ - (i) A x p A ) • + (r 5 - r A ) • (©^ x e A ) = 

2(r s - (r s - x®, xe,) + [(w, xc), x Wj xp^J-e, (12) 

Differentiating once with respect to time in Eq. (8) yields 

(®*xe*)-e^ + e*-(a^xej = 0 (13) 

and differentiating once more gives 

(w s - CO* ) • (e s x e A ) = (w fi - ) • [e^ x (© s x e s ) - e s x x e A )] (14) 

When relative translation is permitted in a certain direction, an equation having the form of Eq. (12) is not re- 
quired. Likewise, when relative rotation is allowed, a relationship in the form of Eq. (14) is not necessary. Instead, 
the required equations can be obtained using the condition that the constraint forces or moments are zero in those 
directions (with the assumption that the surfaces of the joint are perfectly smooth). Thus, for a joint that permits 
translation in a certain direction, 

F(con) . e = 0 (15) 

Depending on the circumstances, F (C0N) can be the constraint force acting either on A or on B. Similarly, when a 
joint permits relative rotation in a certain direction, 


T(con) . e = o 


(16) 


In summary, equations having the form of Eqs. (12) and (14)-(16) provide, in combination, a total of six scalar 
equations. 


The 24 scalar equations formed from Eqs. (l)-(6), (12), and (14)-(16) are linear in 24 unknowns; 3 variables rep- 
resent each mass center acceleration, rigid body angular acceleration, constraint force, and constraint moment for 
each of the two rigid bodies. These equations can be expressed in matrix form, Ax = b (where the column matrix x 
contains the 24 unknown parameters) and solved accordingly using a standard matrix inversion technique. The ma- 
trix A can be singular when, for example, one neglects a central principal moment of inertia, such as for a slender 
rod, or when the mass of a body is neglected. The matrix also becomes singular when constraint equations are re- 
dundant; in other words, when one is linearly dependent on the others. 


III. Test Cases 


A suite of test cases was de- 
veloped to check and verify the 
CFE methodology for application 
to multibody separation problems 
in the launch vehicle staging en- 
vironment. These test cases are 
listed in Table 1 and cover the 
types of joints typically needed to 
model stage separation problems. 
The list is not exhaustive, and, in 
general, it is possible to model 


Table 1. Test Cases Performed for CFE Verification 


Joint Type 

Translational Degrees 
of Freedom 

Rotational Degrees 
of Freedom 

Fixed 

0 

0 

Revolute (Hinge) 

0 

1 

Universal 

0 

2 

Ball 

0 

3 

1-D Translational (Slider) 

1 

0 

2-D Translational (Planar) 

2 

0 

Cylindrical 

1 

1 
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joints that simultaneously constrain any combination of translational and rotational degrees of freedom with CFE 
simply by forming the appropriate equations described in Eqs. (12) and (14)-(16). 

Two test cases from Table 1 (the fixed joint and the universal joint) are described in this paper to illustrate the 
testing and verification that was performed on the CFE methodology. A third test case was conducted that utilized 
POST2/CFE to model the Hyper-X stage separation problem as a 1-D translational (or sliding) joint in series with a 
revolute (or hinge) joint. 

A. Test Case 1: Fixed Joint 

This test case involves two rigid bodies, denoted as bodies A and B , connected by a fixed joint (see Fig. 2). The 
mass properties for each body are shown in Table 2. The fixed joint constrains a point A fixed in A to remain coin- 
cident with a point B fixed in B. Thus, three constraint equations having the form of Eq. (7) can be written as 

(r s -r,)-a,=0 (r= 1,2,3) (17) 

where the role of e A in Eq. (7) is played in turn by , a 2 , and a 3 , each of which belongs to a set of three right- 
handed, mutually orthogonal unit vectors fixed in A as shown in Fig. 2. Suppose that b x , b 2 , and b 3 are a similar set 



Figure 2. Test Case 1: Two rigid bodies, A and B, connected by a 
fixed joint. 


Table 2. Mass Properties for 
Test Case 1, Fixed Joint 



Body A 

Body B 

m 

1.0 kg 

0.4 kg 

Ixx 

0.42 kg-m 2 

0.17 kg-m 2 

Iyy 

2.17 kg-m 2 

0.17 kg-m 2 

Izz 

2.42 kg-m 2 

0.27 kg-m 2 


of unit vectors fixed in B such that b r has the same direction as a r (r = 1, 2, 3) when A and B are attached to each 
other. In that case, the fixed joint constrains certain unit vectors fixed in A to remain perpendicular to other unit vec- 
tors fixed in B. For example, one can write the following three constraint equations having the form of Eq. (8). 


b 2 =0 

(18) 

b 3 =0 

(19) 

b 2 =0 

(20) 


Thus, a total of six constraint equations have been written. As mentioned previously, these constraint equations 
are differentiated twice with respect to time so that the resulting relationships involve the unknown accelerations and 
angular accelerations. Thus, one employs three equations having the form of Eq. (12) in which the role of e A is 
played by a 1? a 2 , and a 3 respectively. Likewise, three equations having the form of Eq. (14) come into play. In the 
first of these, the roles of e A and e B are played by a t and b 2 . In the second such equation, and b 3 play the parts of 
e A and e B . In the third relationship, a 3 and b 2 are substituted for e A and e B . 

No relationships having the form of Eqs. (15) or (16) are applicable in the case of the fixed joint because no rela- 
tive motion is permitted by the joint. 

B. Test Case 2: Universal Joint 

This test case involves two rigid bodies, A and B , connected by a universal joint. This configuration is shown in 
Fig. 3. The universal joint constrains a point B fixed in B to remain coincident with a point A fixed in A, at the inter- 
section of the two arms in the cross. Therefore, the three translational constraint equations presented in Eq. (17) for 
the fixed joint apply in this case as well. 


6 

American Institute of Aeronautics and Astronautics 


For this problem, one additional rotational constraint equation having the form of Eq. (8) is needed. If a 2 is cho- 
sen to be parallel to the arm that is fixed in A, and b 3 is chosen to be parallel to the arm that is fixed in B , then the 
rotational constraint equation imposed by the universal joint can be described with the relationship 

b 3 -a 2 =0 (21) 


The foregoing four constraint equa- 
tions are differentiated twice with re- 
spect to time. Thus, for the transla- 
tional constraint one employs three 
equations having the form of Eq. (12) 
in which the role of e A is played by a t , 
a 2 , and a 3 respectively. Likewise, for 
the rotational constraint the result is 
obtained by employing Eq. (14) and 
making two substitutions: a 2 for e A , 
and b 3 for e B . 

The other two equations are ob- 
tained with the aid of Eq. (16): it is 
assumed that the universal joint is ideal 
(frictionless), so that no constraint 
torque is exerted in the directions of 
the arms in the cross. 



Figure 3. Test Case 2: Two rigid bodies, A and B , connected by a 
universal joint. 


T^(CON) . a 2 = 0, T 5 ( CON )-b 3 =0 


( 22 ) 


The mass properties of each vehicle are summarized in Table 3 
(products of inertia are zero for each body in this problem). 


Table 3. Mass properties for Test 
Case 2, Universal Joint 



Body A 

Body B 

m 

622.0 slugs 

622.0 slugs 

I XX 

220.3 slug-ft 2 

881.2 slug-ft 2 

Iyy 

3524.7 slug-ft 2 

2695.3 slug-ft 2 

Izz 

3330.3 slug-ft 2 

1917.8 slug-ft 2 


C. Test Case 3: Hyper X 

This test case was selected to evaluate the capability of POST2/CFE 
to model and simulate a realistic stage separation problem by compar- 
ing results with flight data from NASA’s X-43A hypersonic scramjet 
test vehicle (see Fig. 4). In 2004, NASA conducted two successful 
scramjet test flights of the X-43A research vehicle (RV) at speeds near Mach 7 and 10, respectively. The roughly 10 
second scramjet test was performed after the RV was boosted to the target flight condition by the Hyper-X Launch 
Vehicle (HXLV). A dimensioned drawing of the mated RV and HXLV is shown in Fig. 5. Additional information 
describing the Hyper-X RV and HXLV may be found in Reference 19. 

The X-43A stage separation event was initiated approximately three seconds after HXLV burnout at a dynamic 
pressure of 1000 psf. Separation began when two pyrotechnically actuated pistons extending from the booster 
pushed against the RV to induce ~16 ft/s of relative velocity between the two vehicles. The two pistons, which were 
positioned roughly 9 in. on each side of the RV centerline, were initially in contact with a cup-like ball joint attached 
to the RV that permitted rotation about the 
piston contact point and translation only 
along the line-of-action of the piston force, 
which was oriented 4 deg below the hori- 
zontal. Both pistons were connected to the 
same gas chamber and were activated at the 
same time and remained in contact with the 
RV for ~0.1 sec until they reached the end 
of their 9-in. stroke length. Although the ball 
joint permitted rotation in any direction at 

each individual piston contact point, the fact separation from HXLV booster. 
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that both pistons were fired concurrently and were simul- 
taneously in contact with the RV effectively excluded 
relative rolling and yawing rotation between the two sepa- 
rating vehicles. 

Once the pistons reached the end of their 9-inch ex- 
tensions, the RV separated from the HXLV. At this time, 
the HXLV commanded a nose-down pitching moment to 
move away from the RV, and the RV was steered to its 
required target condition. 

To test the CFE methodology, this problem was set up 
using POST2/CFE to model the piston constraint forces. 
For the purpose of this paper, the two pistons were mod- 
eled as a single piston to avoid over-constraining the 
problem. Thus, a single piston was modeled as a sliding 
joint in series with a revolute joint; that is, relative trans- 
lation was allowed in one direction (the piston line-of- 
action) and relative rotation was permitted only about the 
pitch axis. The piston was placed at the centerline of the 
RV and the axial force that pushed the vehicles apart was 
doubled. A schematic highlighting these modeling details 
is shown in Fig. 6. 

To illustrate how the constraint equations for this joint 
are derived, unit vectors are defined in the RV such that 
a 1 is the direction in which translation or sliding is permit- 
ted (along the piston line-of-action) and a 2 indicates the 
direction in which rotation is free to occur (pitch axis). A 
similar set of unit vectors, b r , can be defined that are 
fixed in the HXLV and initially aligned with the a r unit 
vectors fixed in the RV (see Fig. 6). For this case, two 
constraint equations having the form of Eq. (7) can be 
written as: 



Figure 5. Test Case 3: Dimensioned drawing of 
mated X-43A Research Vehicle and Hyper-X 
Launch Vehicle. 




Figure 6. Test Case 3. CFE modeling of piston 
contact for Hyper-X (X-43A) separation problem. 


(r fi -r,)-a r =0 (r = 2,3)) (23) 

As before, these two constraint equations are differentiated twice with respect to time as in Eq. (12). Now, s^and 
a 2 are substituted for e A . Furthermore, two constraint equations describing the restriction imposed on relative orien- 
tation are formed, 


a 2 • b, = 0, a 2 • b 3 = 0 


(24) 


Similarly, these constraint equations are differentiated twice with respect to time as in Eq. (14), where now a 2 is 
substituted for e A , andbj and b 3 are substituted for e A . With the assumption that the joint is frictionless, the con- 
straint force and moment in the directions where motion is permitted become, using Eqs. (15) and (16), 


|7(C0N) 


•a, =0, 


t (CON) 

l B 


a 2 = 0 


( 25 ) 
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IV. Results and Discussion 


A. Test Case I, Fixed Joint 

The two rigid bodies, A and B , are assumed to be rigidly connected to each other for 10 sec (see Fig. 2) and then 
released instantaneously. No external forces or moments were assumed to act on either body for the entire duration 
of the 20-sec simulation. Initial conditions correspond to zero inertial velocity of the system mass center (£*), and an 
inertial angular velocity of the rigid body formed by A and B equal to 63.02^ -82.32a 2 +80.25a 3 deg/s. After 10 
sec, the joint connecting A and B is released, and their motions become unconstrained. Because there are no external 
forces or moments, integrals of the motion involving linear and angular momentum must remain constant and this 
fact can be used to check numerical solutions of the equations of motion. Therefore, one can expect the simulation 
results to show that, after separation, the mass centers of A and B travel in straight lines with constant velocities, 
whereas the system mass center remains at rest. 

To verify the POST2/CFE simulation results, a simulation was also created with AUTOLEV, 20 " 22 an interactive 
program designed specifically for the kinematic and dynamic analysis of mechanical systems. Symbolic manipula- 
tion (computer algebra) is used to formulate explicit equations of motion for the particular system of interest using 
any method; the software is particularly well suited to Kane’s method. 22 ' 23 For this test case, AUTOLEV was used to 
create a computer program to simulate the motions of A and B , and determine the constraint forces required to hold 
them together for the first 10 seconds. 

The AUTOLEV results were generated using a variable step integrator with an absolute error limit of 1 x 10” 8 
and a relative error limit of 1 x 10” 7 . The POST2/CFE method employed a fixed step integrator with a step size of 
0.0001 sec. The CFE routine also used a Baumgarte 24 constraint stabilization factor of 5.0. 

Time histories for the inertial x, y, and z components of velocity of the mass center of body A are shown for 20 
sec in Fig. 7. The values oscillate during the first 10 sec, after which the joint is released and their values become 
constant, as expected. Figure 7 shows that there is an excellent agreement between POST2/CFE and AUTOLEV 
results. Similarly, the velocity components for the mass center of body B in Fig. 8 also indicate an excellent agree- 
ment between POST2/CFE and AUTOLEV. 

The angular velocity time history comparisons for bodies A and B are shown in Figs. 9 and 10, respectively. 
Again, there is excellent agreement between the two approaches. The change in the angular velocity of body B at the 
point of release is more noticeable than for body A because B is smaller and has a mass distribution that differs more 
from that of the composite body before joint release. The main results of the CFE method are the constraint forces 
and moments, which are shown in Figs. 11 and 12, and are applied as additional external loads within POST2. The 
figures show that the forces and torques determined by CFE agree very well with the AUTOLEV results. 

An important metric to assess the POST2/CFE methodology is the relative joint displacement between the two 
bodies when they are supposed to stay connected. This parameter is computed as the position vector from A to B and 
should be zero while the joint constraint is imposed for the first 10 sec. Hence, any deviation from zero is a measure 



Figure 7. Test Case 1, Body A: Velocity components 
of mass center. Joint release occurs at t = 10 sec. 


POST2/CFE as circles, AUTOLEV as lines. 



Figure 8. Test Case 1, Body B: Velocity components 
of mass center. Joint release occurs at t = 10 sec. 
POST2/CFE as circles, AUTOLEV as lines. 
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components. Joint release occurs at t = 10 sec. 
POST2/CFE as circles, AUTOLEV as lines. 



Figure 11. Test Case 1: Comparison of constraint 
forces required to hold bodies together. Joint release 
occurs at t = 10 sec. POST2/CFE as circles, 
AUTOLEV as lines. 



components. Joint release occurs at t = 10 sec. 
POST2/CFE as circles, AUTOLEV as lines. 



Figure 12. Test Case 1: Comparison of constraint 
torques required to hold bodies together. Joint release 
occurs at t = 10 sec. POST2/CFE as circles, 
AUTOLEV as lines. 


of accuracy of the CFE algorithm. In the CFE algorithm the constraints are satisfied identically at the acceleration 
level [Eqs. (12) and (14)]; however, as with all algorithms of this type, the constraints at the velocity and position 
levels are inevitably subject to numerical integration errors. Such errors are functions of the step size in a fixed-step 
integration scheme or the error limits in a variable-step approach. Baumgarte stabilization 24 is used to control these 
errors from growing arbitrarily large during a simulation. 

The choice of step size and Baumgarte factor for this problem represents a reasonable balance between CPU 
time and error buildup. Figure 13 shows the joint translational displacement as a function of time for the 10 sec pe- 
riod when the bodies are supposed to stay connected. In this paper, the quantities Ax, Ay, and Az are the x, y and z 
components of joint translational displacement in the local body frame and are defined as follows: 

Ax = (r s -r j4 )-a 1 =0, Ay = (r s -r„)-a 2 =0, Az = (r B -r^)-a 3 = 0 (26) 


In the interest of clarity, only the POST2/CFE results are shown. The magnitude of this displacement distance is 
below 0.01% of the total body length. 


10 

American Institute of Aeronautics and Astronautics 



Figure 13. Test Case 1: Joint displacement in x, y 
and z directions for POST2/CFE method during 
time when bodies are connected. Joint release 
occurs at t = 10 sec. 



Figure 14. Test Case 2: External force components 
applied to Body A (x and y components are identical). 


B. Test Case 2: Universal Joint 

To check and verify POST2/CFE results, an independent simulation using ADAMS 12 software was generated. At 
t = 0, each body was assumed to be at rest. A time-varying force is applied to Body A at a location relative to the 
center of mass given by the vector -2.0^ + 0.5a 3 . The components of this time-varying force are shown in Fig. 14. 
The mass properties of each vehicle are summarized in Table 3. The products of inertia are zero for each body in 
this problem. 

The motion of bodies A and B was simulated for a duration of 5 sec. The POST2/CFE used a fixed step size in- 
tegrator with step size of 0.0001 sec with a Baumgarte factor of 5.0. The results of the simulation cases are shown in 
Figs. 15-18. The angular velocity of body A is shown in Fig. 15. The constraint force and moment components are 
shown in Figs. 16 and 17, respectively. There is excellent agreement between POST2/CFE and the ADAMS solu- 
tions as evidenced by these figures. Finally, the relative joint translational displacement time histories are shown in 
Fig. 18. As in the previous test case of the fixed joint problem, the joint displacement is relatively small. 



Components. POST2/CFE as circles, ADAMS as constraint force. POST2/CFE as circles, ADAMS as 
lines. lines. 
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Figure 17. Test Case 2, Body B: Comparison of 
constraint torque. POST2/CFE as circles, ADAMS as 
lines. 


Figure 18. Test Case 2: Joint displacement in x, y 
and z directions for POST2/CFE method during 
time when bodies are connected. 


C. Test Case 3: Hyper X 

The results of POST2/CFE simulations were compared to previous POST2 simulations based on an engineering 
joint model using springs 25 and the flight data. 26 The engineering model employed massless springs to model sliding 
joints at both piston locations. As a measure of how well the joint constraint was satisfied, the relative joint dis- 
placement between the two separating vehicles is shown in Fig. 19. This parameter is the distance from the RV joint 
location to the HXLV joint location, and was computed in the piston axial (line-of-action, aj, lateral (a 2 ), and 
normal ( a 3 ) directions. The results indicate that CFE does a better job of enforcing the translational constraint, and 
there is virtually no joint displacement in the lateral or normal directions. Similar results are seen in Fig. 20 which 
shows the relative difference in attitude between the two vehicles. Again, CFE does an excellent job of enforcing the 
rotational constraint, essentially eliminating any relative roll and yaw between the RV and HXLV. The constraint 
forces and moments about the vehicle center- of- gravity that were computed by each technique are shown in Figs. 21 
and 22, respectively. In these plots, the forces and moments computed by the spring model include the contributions 
of both pistons. The results show that both the spring model and POST2/CFE compute very similar forces and mo- 
ments. However, the spring model results are oscillatory in nature since the joint constraint is not satisfied exactly. 


Spring Model 

CFE 

o RV Free (end of piston extension) 



0002 1 . 

d 0.001 - I A - 

< -o.ooi - i i 

-0.002 i i 

0 0.05 0.10 0.15 

Time, sec 



Figure 19. Test Case 3: Comparison of joint dis- 
placement in x, y and z directions CFE and Spring 
Model. 



Time, sec 

Figure 20. Test Case 3: Comparison of relative 
angular displacement between POST2/CFE and 
Spring Model. 
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Figure 21. Test Case 3: Comparison of constraint 
forces computed by POST2/CFE and Spring Model. 



Figure 22. Test Case 3: Comparison of constraint 
moments computed by POST2/CFE and Spring 
Model. 
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sideslip profiles are shown in Fig. 23. Both the spring — cfe, adj. aero 

11 1 i i 11 i Unconstrained, adj. aero 

model and POST2/CFE results match very well and are o Flight Data 

nearly identical in the angle-of-attack response. There are 
slightly larger differences in the sideslip response, but most 
of the motion occurs in the pitch plane. In addition, both 
simulation techniques match the observed flight data well. 

It should be noted that the degree to which the simulation 
results match the flight data depends on other factors be- 
sides CFE modeling, such as aerodynamic force and mo- 
ments, winds, flight data instrumentation, etc. Aerody- 
namic modeling was especially critical because separation 
occurred at a dynamic pressure of 1000 psf. Thus, this 
problem is much more complex than the first two test prob- 
lems where the external forces and moments were known 
exactly. 

To illustrate the sensitivity of the simulation results to 
aerodynamic modeling, results are included using the 
nominal (pre-flight) aerodynamic database and a post-flight 
reconstructed aerodynamic database in which the uncer- 
tainty terms had been adjusted to better match the flight 

data. The variations in angle-of-attack and sideslip differ significantly depending on which aerodynamic model was 
used. Note that all POST2/CFE results were generated using the adjusted aerodynamics. 

As a further comparison, results are also shown for an unconstrained case where no constraint forces were ap- 
plied to model the joint (but the piston axial force that pushed the vehicles apart was still modeled). During the first 
second of separation, the unconstrained results differed significantly from the cases that included the constraint, 
demonstrating the importance of accurate modeling of the piston joint. Even though the differences between the 
constrained and unconstrained angle-of-attack and sideslip profiles disappeared after the control system became 
fully active and was able to compensate, it is important to capture the dynamics during the first second of flight 
since that was when the risk of re-contact was highest. 

Similar comparisons are shown for the linear accelerations in Fig. 24. The vertical (Body Z direction) accelera- 
tion results were the most sensitive to the way in which the piston constraint was modeled, and was also the most 
important for assessing the risk of re-contact. Again, the spring model and CFE results were nearly identical and 
provided a good match to the flight data. The differences between the constrained and unconstrained cases were 
actually larger than the differences due to aerodynamics. 


Figure 23. Test Case 3: Simulation vs flight data 
comparison of angle-of-attack and angle-of-sideslip 
profiles. 
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Figure 24. Test Case 3: Simulation vs flight data 
comparison of x,y and z accelerations in the local 
body frame. 




comparison of body roll, pitch and yaw rate 
profiles. 


Figure 25 compares the angular velocity results with the flight data. Once again there is excellent agreement be- 
tween the spring model and CFE simulation results. Both simulations match the pitch rate flight data very well and 
do a much better job at predicting the pitching motion than the unconstrained case. The comparison between simula- 
tion and flight for yaw rate and roll rate is not as good. Much of the difference can be attributed to inadequate lat- 
eral-directional aerodynamic modeling as is evident in the large difference in results between the cases with nominal 
aerodynamics and adjusted aerodynamics. Again, capturing the pitch plane dynamics was important for assessing 
the risk of re-contact. 


V. Concluding Remarks 

A suite of test cases was designed to verify the POST2/CFE methodology and to assess its suitability to simulate 
multibody separation problems in the launch vehicle staging environment. To illustrate the capability of the CFE 
algorithm that was implemented in POST2, three test cases of increasing complexity were presented. The first case 
examined the motion of a two unperturbed rigid bodies connected by a fixed joint subject to an initial rotation. The 
simple nature of this problem made it possible to compare simulation results to analytical predictions based on first 
principles. In addition, simulating a fixed joint presented a challenging case in which to assess the relative displace- 
ment between the two bodies at the joint location, a key performance metric for the CFE algorithm. The second case 
involved rigid bodies connected by a universal joint and typified the testing that was performed on a number of dif- 
ferent joints applicable to stage separation problems. For the first test case, the POST2/CFE results were in excellent 
agreement with results obtained using AUTOLEV and for the second case, with those obtained using ADAMS. 

In the third test case, the ability of POST2/CFE to model a realistic launch vehicle stage separation problem was 
demonstrated by simulating the X43-A (Hyper-X) stage separation problem. The POST2/CFE results agreed well 
with previously published POST2 results that were based on an engineering model and actually performed better at 
satisfying the joint constraint. Moreover, the POST2/CFE results agreed reasonably well with the flight data, and it 
was shown that the level of agreement degraded significantly if the constraint was not modeled. 

The success of these test cases provides confidence in the ability of the CFE methodology to solve constrained 
motion problems in conventional trajectory simulation programs. The CFE implementation in POST2 provides the 
capability of simulating generic end-to-end simulations of launch vehicle trajectories. 
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